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The concept of transmission eigenchannels is described in a tight-binding nonequilibrium Green's 
function (NEGF) framework. A simple procedure for calculating the eigenchannels is derived using 
only the properties of the device subspace and quantities normally available in a NEGF calculation. 
The method is exemplified by visualization in real-space of the eigenchannels for three different 
molecular and atomic-wires. 
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I. INTRODUCTION 



Electronic transport properties of atomic-scale 
conductors has been investigated intensively in the 
last decadeii^. Examples of interest include molec- 
ular wires connected to metal electrodes, atomic 
metal wires, and nanotubes. First principles trans- 
port calculations on these systems give results that 
are, in general, difficult to interpret due to the multi- 
channel nature of the scattering problem and the 
fact that the scattering states are generated from 
the atomic valence orbitals. The free-electron-type 
of reasoning normally used in mesoscopic quantum 
transport is thus not adequate. It is therefore use- 
ful to analyze the conduction in terms of eigen- 
channels. Eigenchannels are particular scattering 
states^ with a well-defined transmission probabil- 
ity, < T„ < 1, where the individual eigenchan- 
nel transmissions add up to the total transmission 
T = X^n-^"- addition to being useful for analyz- 
ing theoretical calculations, the eigenchannel trans- 
missions may be obtained experimentally (i) with 
superconducting electrodes connecting the atomic- 
scale conductor, as first shown by Scheer et ali^, or 
(ii) from shot noise measurements^ where informa- 
tion about the individual channel contributions can 
be obtained since the Fano factor involve the sum^'^ 

E„T„(i-r„)/E„T„. 

Eigenchannels have previously been calculated for 
atomic metal wires^ and molecular contacts^ by di- 
rectly solving for the scattering states in the leads. 
This type of analysis breaks up the transmission into 
the 'non-mixing' channels^i^iii^ and gives an intuitive 
picture of electron transport. The ability to plot the 
eigenchannel wavefunctions is especially useful, since 
it gives a direct spatially resolved picture of the or- 
bitals involved in the transport. Another possibility 
is to consider projections of eigenchannel wavefunc- 



tions onto, for example, molecular orbitals. 

An increasingly popular theoretical approach to 
calculate transport properties is the non-equilibrium 
Green's function (NEGF) formalism^^. This is nor- 
mally used in combination with a tight-binding type 
or LCAO electronic structure description— i^'^i^^i^^ . 
In this approach, it is straightforward to calculate 
the single particle (retarded) Green's function (ma- 
trix) including coupling to the infinite electrodes by 
introducing self-energies. The Green's function is 
thus the fundamental quantity in these calculations 
and scattering states are normally not considered. 
Interpretation of the results in terms of scattering 
states is therefore non-triviaU^. In contrast, the 
scattering states are the fundamental quantity in 
approaches based on the Lippmann-Schwinger equa- 
tion where jellium models are normally used to de- 
scribe the electrodesiii^ii^'SS. 

The aim of this paper is to show how the eigen- 
channels can be easily generated within the NEGF 
approach without solving for the scattering states 
in the leads. The eigenchannel wavefunctions are 
here obtained directly from quantities readily avail- 
able in the NEGF calculation, e.g., the retarded 
Green's function matrix, G^i of the device region, 
and the Tlh matrices describing the coupling of 
the device region to the two electrodes ('left' L and 
'right' R, see Fig. [T]). In the case of atomistically 
defined electrodes, this approach is especially ad- 
vantageous since solving for the scattering states 
requires calculating the Bloch waves in the elec- 
trodes (complex band structure), which may be a 
non-trivial numerical task for large unitcells^. Re- 
lated to our approach is the so-called (left/right) 
open 'channel functions' of Inglesfield an co-workers 
based on the 'embedding potential' in the real-space 
formulationSS,, and the KKR-based formulation by 
Bagrets et ali^^^. In addition to providing a sim- 
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pie way to calculate the eigenchannels, the method 
presented here offers an intuitive understanding of 
the one-particle NEGF equations, e.g., it may be 
used to understand propensity rules for the effect of 
phonon scattering on the electronic transport^^'^^. 

The paper is organized as follows. Sec. [Til starts 
with the definition of eigenchannels and a summary 
of the standard one particle NEGF equations. Our 
method to obtain the eigenchannels without solv- 
ing for the full transmission matrix is then derived. 
The usefulness of the eigenchannels is illustrated in 
Sec. mil with three examples where the eigenchan- 
nels are calculated for molecular and atomic-wires 
connected to gold electrodes using a first-principles 
density functional method. 

II. EIGENCHANNELS 
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FIG. 1: The generic two probe system which couples 
the left (L) and right (i?) leads through an intermediate 
"device" region (D). 

We consider transport through a device region, 
Z?, coupled to two semi-infinite leads, left and right 
(L, R) and limit ourselves to leads built from peri- 
odic cells. The solutions in the corresponding infi- 
nite left /right leads are Bloch states \ui) (|ur)) de- 
noted by band-index I (r) for the left (right) lead. 
In order to obtain the transmission amplitude ma- 
trix t,..; at a given energy, E, we consider the solu- 
tions to the Schrodinger equation, |^'/) with scatter- 
ing boundary conditions and incoming waves in the 
left lead, 

Y^i) = + Reflected ; left lead, 

1*0 = E 7=tr,; + Decaying ; right lead, (1) 

where \ui) is an incoming Bloch wave from the left 
lead, i.e., right-moving, with energy E and group 
velocity vi > 0, and \ur) is an out-going Bloch wave 
{vr > 0). At a certain energy, the number of such in- 
coming channels, Nl(E), is determined by the band 
structure of the lead. Likewise, there are Nji(E) 
such out-going channels on the right. In Eq. [U we 



explicitly state the flux-normalization by dividing 
the Bloch state by y^, where the Bloch waves \ui) 
are normalized in the conventional manner over the 
infinite leadsSi {ui{k)\ui>ik')) = 5w 5{k - k') . The 
reflected and transmitted parts may contain evanes- 
cent decaying waves which have zero velocity and 
for those states we use normal integral normaliza- 
tion. With these considerations we find that the 
flux-normalized states fulfil the normalization, 

-^{ui{k)\ui,{k'))^^ hSn'S{E~E'), (2) 

since the velocity is related to the energy as vi = 
l/hdei{k)/dk. The scattering states generated from 
the flux-normalized incoming waves will also be flux- 
normalized, see Appendix [XI 

{^i\^i,)^h6ivS{E~E'), (3) 

and likewise for r, r' while = 0. 

The advantage of using flux-normalized states is 
that we can make any unitary transformation be- 
tween the incoming scattering states from the left 
lead at a particular energy, 

K)=5:^(u^).,. (4) 

These new states will again solve the Schrodinger 
equation, only contain incoming waves in the left 
lead, and be fluxnormalized. However, the mix of 
Bloch waves will no longer have Bloch symmetry. 
Naturally we can apply a similar transformation of 
the out-going right channels with Vr. Especially 
we can choose the transformations (Ul,Uk) such 
that the transmission matrix t becomes a diagonal 
matrix^, te, at a specific energy 

t, = Uj^ t Ui - Diag (v/l\, Vt'2, ■■)■ (5) 

This corresponds to a singular value decomposition 
(SVD) of the transmission amplitude matrix, or a di- 
agonalization of the Hermitian (left-to-right) trans- 
mission probability matrix, 

Te =ttte = u[tttUL =Diag(Ti,r2,...) , (6) 

where there will be a maximum of 
Mm{NL{E), Nii{E)) non-zero, eigenvalues 

< Ta < 1- 

The diagonalization of the transmission matrix 
defines the transmission eigenchannels \^a) of the 
system as the unitary mix of left in-coming flux- 
normalized scattering states ("channels") given by 
Ul, i.e., the left eigenchannels are given by \^a) ~ 
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J2 |^'/')(UL);',a. The eigenchannels have the special H\'i'.m{E)) = E\'i'„^{E)), 



property of being "non-mixing" in the sense that the 
transmission of a sum of two is the sum of individ- 
ual fluxes or transmissions (equal for flux normalized 
states). For example, consider the alternative chan- 
nels given by the scattering states defined by the two 
first eigenchannels as, \ua) = ai|$i) +a2|<l'2), with 
transmission Ta — |aipri + |a2p72 and similarly \ui,) 
with transmission Tf,. Now the total transmission of 
\uc) — \ua) + \ub) will be Tc — Ta+Tt+I, with the in- 
terference term / = 2TiRe(ai6i) + 2 T2Re{a2b2) ■ It 
will thus not simply be the sum of the two transmis- 
sions, Tc y^Ta + Tf,. Since the purpose of this paper 
is to calculate the eigenchannels without solving for 
the complex band structure in the leads, we will use 
the fact that the eigenchannels also maximize the 
transmission through the device, i.e., the first eigen- 
channel from the left contact maximize the trans- 
mission probability over the space of incoming states 
from the left, the next channel maximizes the trans- 
mission while being orthogonal to the first channel 
etc. 

In the following, we will instead of flux- 
normalization make use of energy-normalization 
(*„(T;) I *„(£;')) = SnmS{E - E') with the trivial 
difference from flux-normalization being a factor Vh. 
Energy-normalization is advantageous when work- 
ing with energy resolved quantities since the natu- 
ral continuous quantum number for this normaliza- 
tion is the energy. In addition, we can interpret the 
energy-normalized states as a density-of-states, i.e., 
|*;(2:)P — is the projected (local) density- 

of-states at X. 



A. Preliminaries 

The one particle Hamiltonian for the tight-binding 
scattering problem shown in Fig. [T] can be written 

H = Ho + V = Hl+Hd + Hr + Vl + Vr, (7) 

where the isolated leads Hl r and device Hu are 
coupled together by the interactions between leads 
and device {V — Vl + Vr) without any direct cou- 
pling between the leads. Using projection operators 
onto the device Pd and left/right leads Pl,r {I = 
Pl+Pd+Pr) we can also define tl,r = PdVl,rPl.,r 
where Vl,r — tl,r + r- The rest of this section 
provide a short summary of the standard definitions 
of one particle Green's functions, self-energies and 
the notation we will use throughout the paper— i^l. 

From the definition of the retarded Green's func- 
tion (operator) for the whole system we can find the 
expansion of the Green's function in an eigenbasis. 



G{E) = [E + iS-H] 



\^,n{E')){^„AE')\ 

E + i5-E' 



(8) 



where the infinitesimal imaginary part 5 = Q'^ en- 
sures that the Green's function yield the retarded re- 
sponse of the system. The device part of the Green's 
function can further be written 



GD = {E + i5-HD-y^L-^Ry 



(9) 



where we have introduced the self-energies, Sl.h = 
''l r.9l R^l. R given by the Green's functions of the 
isolated leads gL,R — {E + iS — Hl,r)~^- In ad- 
dition to the Green's function, the spectral func- 
tions A{E) = i{G-G^), aLM = i{9LM- sI^r), and 

broadening Fl.^r = i{^L,R - ^l,r) = TL,RaL,RTl R 
will be needed. The fact that these matrices live 
on different subspaces will be used repeatedly, e.g., 
Tl = PdTlPd etc. 

For the scattering problem, see Fig. [U we know 
that the time independent discrete Schrodinger 
equation has a complete set of solutions. These so- 
lutions can be divided into a continuous set of solu- 
tions \'^n{E)) (where there may be several solutions 
at any given energy) and localized states 1^'^^°'^) with 
energy We use the energy as the continuous 

quantum number together with a discrete quantum 
number n, i.e., sub-bands. Since we are only in- 
terested in the transport properties, we will from 
hereon ignore localized states^. 

In the following, we will describe a method to 
determine the eigenchannel scattering states inside 
the device region using the information contained 
in G_D, Fl, and Tr. The spectral function, A, is 
a central quantity in the following discussion. It 
can be obtained from the expansion of the retarded 
Green's function in eigenfunctions to the Hamilto- 
nian (Eq.[5D, 

AiE) = ^{GiE)~GHE)) 

= 2nJ2\'^n{E)){^n{E)\. (10) 



B. Scattering states from the leads 

We may choose to express the solutions to the 
Schrodinger equation as solutions consisting of waves 
originating in the left or right lead. These scattering 
states can be generated from the spectral function as 
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will be show here. Decomposing the spectral func- 
tion of the device Ad — PdAPd, using Eq. ([9]), we 
find 

Ad = i{Gd - G^,) - iGd (g^d^^ - G^^) G^, 

— GuF^G^ + Gd^rG^jj , (11) 
where we in the following wish to show that Aj^ n = 

Gd^l,rG]j is generated by the scattering states 
with incoming waves in the left (right) lead. 

Viewing the coupling between device and leads as 
a perturbation, we can start with a set of orthogonal 
and normalized eigenfunctions, \ui), of the isolated 
left lead (and similarly for the right) which are to- 
tally reflected solution since they are solutions for 
the isolated semi-infinite leads. From these states, 
the full solutions can be generated, 

=GVl\ui) + \ui) . (12) 

The response given by the retarded Green's function 
only contains waves traveling outwards from the de- 
vice region. This solution to the Schrodingcr equa- 
tion thus have the required property of being incom- 
ing from the left lead. In addition, the solutions are 
energy-normalized and orthogonal, see Appendix \X[ 
We can then express the device part of the spectral 
function from the solutions generated by Eq. (fT2|) 

Al - 27rJ2 PD\'fi){'i'i\PD 

= 27rJ2 Pd{GVl\ui) + \ui)) 
I 

[{ui\ + {ui\vIg^)Pd 

= PdGVlulVIg^ Pd = G^FiG^, , (13) 

where we have used Eq. (fTO| for the whole system 
and for the isolated lead. Apart from re-deriving 
Eq. (fTTj) we have proven that the two parts of the de- 
vice spectral function Aj^^^i — Gd^l.iiG\j are built 
up of scattering states originating from the respec- 
tive leads. This immediately leads to the well known 
formula for the density matrix in non-equilibrium 
(excluding localized states) 

oo 

P=^J [fLiE)AL{E) + fRiE)AiiiE)] dE, 

— OO 

(14) 

where /l,_r is the Fermi function of the leads. 

C. Current operator 

To find the eigenchannels of the system we will 
need the current operator. The number of electrons 



in lead R is described by the projection operator Pr. 
The operator for current into R is thus determined 
the time derivative of Pr, 

j^ = 2ePn='^[H,Pn] = '^{r„-rl,), (15) 

where we have evaluated the commutator using the 
Hamiltonian, Eq. ([7]), and included a factor of 2 for 
spin. 

The current into lead R due to the scattering state 
with energy E, \'i>i), originating from L from the 
original in-coming (standing) wave \ui), can be writ- 
ten 

M^{^i\J\^i). (16) 

To simplify this expression, we extract the right 
lead part of the wavefunction using the Lippmann- 
Schwinger Eq. (|A2[) . gives 

PrI^i) = Pr.GoV\^i) + Pr,\ui) = cjrtII^i) , (17) 

where all quantities are evaluated at energy E. The 
current carried by the scattering state is then 
given by 

= ^ 2n{^i\rR\-^i} . (18) 

Summing the current over all the orthogonal, 
energy-normalized scattering states originating in L 
at the specific energy yield the net current or total 
transmission, and thus the Landauer formula, 

J = E'^" = ^27r^(vl/,|rfl|vl/,) (19) 

= Tr (AlTr) = ^ Tr {GdTlG^oTr^ . 

From Eq. [THl we notice that the transmission 
probability for any scattering state from the left \^i) 
is given by 2tt {'i> i\T r]'^ i) . The transmission proba- 
bility matrix can therefore be written 

Tin ^27r{^ I, ITrI^i) , (20) 

since we consider energy-normalized (flux- 
normalized except for the factor of h) scattering 
states, l^*), for which the current is equivalent with 
the transmission. 
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FIG. 2: Overview of the different states used in this sec- 
tion. Note that the scattering states which are non-zero 
in the device subspace are aU spanned by the |^';) states. 
We can therefore find the eigenchannels by studying the 
wavefunctions spanned by the |xi) = PdY^i)- 



D. Scattering states in the device region 



To find the eigenchannels from the left lead, |<i>;) = 
1^'/'), we need to diagonalize the transmis- 

sion probability matrix. Since we do not, in this pa- 
per, explicitly calculate the Bloch states in the leads, 
we will diagonalize the transmission probability ma- 
trix in an abstract basis formed by incoming waves 
from the left or right lead. Another, equivalent, for- 
mulation is to find the eigenchannels by maximizing 
the transmission at a given energy while keeping the 



eigenchannels orthogonal (UiU^ = 1) 
max 27r(<i>;|r7? |<i>;) . 



(21) 



The main problem in finding the eigenchannels is 
that we do not have access to the wavefunctions. 
Green's functions or spectral functions of the entire 
system in a typical calculation. We therefore have to 
find a method using only properties from the device 
part of the system. 

We will now show that optimizing the current by 
varying the scattering state over space spanned by 
the incoming states from the left lead, {1x1/;)}, is 
equivalent with varying over the states \xi) defined 
on the device subspace, see Fig. [51 We start by di- 
agonalizing the device part of the spectral function 
from the left lead 



Ai(i?) = ^|Xi)A,tol = 2^5^|Xi)(x/ 



(22) 



where the eigenvectors on the finite device space 
are orthonormal {xi\xi') = 5i,i' ■ Each \xi) = 
^J\i/2'k\xi) with non-zero eigenvalues (A; 7^ 0) is 
the device part of a specific state |^;), i.e., \xi) = 
^bl^i), where the |^';) states are normalized and or- 
thogonal linear combinations of the scattering states 



\xi) 




■■Al{E)\xi) 



(using Eq. (fT3|) ) where we have defined 
1^') = E a/|^(*/'Ix;)I*p) ^ E I*'') ■ (24) 



This shows that the states \'^i) are spanned by the 
incoming scattering states from the left lead. We 
further wish to show that these states are normalized 
in the same manner as the original scattering states, 
i.e., we want to show that Wmn is unitary 
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(25) 



using Eqs. ^ and (j^ . 

In addition, any scattering state outside the space 
spanned by \^n) are orthogonal to the device sub- 
space. This can be seen from the fact that we can 
write the spectral function on the device subspace 
as 

Al = 2t:^Pd\^i)1^i\Pd or as (26) 

I 

= 27: Y,PD\'fi){'i'i\PD- (27) 
I 

Comparing the two equations reveals that 
P]j\'^){'i'\P]:) must be zero for any j^*) which 
is orthogonal to the space spanned by l^*;). 

In this section we have shown that the wavefunc- 
tions Ixi) span the device part of any scattering 
state generated from lead L, see Fig. To find the 
eigenchannels we can therefore maximize the current 
through the device with respect to a linear combi- 
nation of \xi) instead of maximizing with respect to 
the full scattering states 1 5';) . This is equivalent with 
diagonalizing the transmission matrix (Eq. (|20[) ) in 
the basis formed by \xi)- 

E. Finding the eigenchannels 

In general, the basis {le^}} used in a calculations 
is non-orthogonal with the overlap matrix defined by 
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Although the eigenchannels may be 



calculated directly in this non-orthogonal basis, we 
will make use of Lovdin ortogonalization to simplify 
the algebra. The ortogonalized matrices (denoted 
by a bar) are given by, Tr = S~^/^ri^S~^/^ and 
Al = S-'^/^AlS^/^ etc. The eigenchannels obtained 
in the Lovdin orthogonalized basis can at the end 
of the calculation simply be transformed back into 
the non-orthogonal basis for further visualization or 
projection. 

From the previous section we learned that it is 
enough to diagonalize the transmission probability 
matrix (Eq. = 27r(x;' |f ijjxi) in the ab- 

stract basis {|xi)}- To transform 2TTTii into the basis 
{|Xi)}i we first need to calculated the eigenvectors 
of AL=GfLG^, 



ii ' 



(28) 



where U is unitary. We then obtain the transforma- 
tion matrix to the {|x;)} basis. 



U 



7nl 




ml ' 



(29) 



which gives the explicit expression for the matrix we 
want to diagonalize 



Tiq = 2'k{xi'\Tr\xi) = 27r 
the eigenproblem is therefore 



I'l 



^27r [Ut] 



(30) 



(31) 



where the eigenchannel vectors, Cq, are given in the 
basis described by the columns of U, and the eigen- 
values Ta are the transmission probabilities of the 
individual eigenchannels, a. Finally, transforming 
back to the original non-orthogonal basis (from U- 
basis to the Lowdin-basis, and from Lowdin basis to 
normal non-orthogonal basis) we find that the eigen- 
channels on the device subspace are given by, 



shown using the formalism presented above). This 
provides a simple method to obtain an idea about 
the eigenchannels. However, it is important to real- 
ize that the eigenchannel wave-functions calculated 
from Eq. [32] are energy normalized, i.e., the ampli- 
tudes are well defined and can be compared between 
different eigenchannels. In contrast, the eigenvec- 
tors to GTiG^Fi? may have any normalization and 
it is therefore not possible to compare amplitudes 
between different channels. Moreover, the energy 
normalized scattering states Eq. [32] yield amplitudes 
which correspond to local density of states and are 
useful to plot, as will be shown in the next section. 



III. EIGENCHANNELS FOR ATOMIC- AND 
MOLECULAR WIRES 

To exemplify the method developend in the pre- 
vious section we will use three different examples of 
atomic and molecular-wires connected to gold elec- 
trodes. Using the TranSIESTA^^ extension of the 
SIESTAS^ density functional theory (DFT) code, 
we have previously studied elastic and inelastic 
transport properties for the systems under consid- 
eration here, (i) atomic gold wirea^^, the conju- 
gated organic molecules (ii) oligo-phenylene viny- 
lene (OPV)^, and (iii) oligo-phenylene ethynylene 
(OPE)30. The TranSIESTA calculations on these 
systems were performed using DFT in the gener- 
alized gradient approximation using the PBE func- 
tional. The semi-infinite leads connecting the device 
region was modeled using self-energies in the NEGF 
method. A more detailed description of the calcula- 
tional method may be found in Ref.— . 




PD\<^a} - 



7n.n 



(32) 



The Eqs. (|28| - (|32]) provide a recipy for calculating 
the eigenchannels of a specific scattering problem us- 
ing only properties available in standard NEGF cal- 
culations. In contrast to Ref.— ^ these eigenchannels 
are well defined scattering states calculated without 
approximations on the full device subspace. 

It is interesting to note that the eigenchannels 
(Eq. [32]) are eigenvectors to GTlG^Tr (can be 



FIG. 3: (Color online) Left eigenchannel for a 7-atom 
An atomic chain connected by 4 atom pyramids to An 
(100) surfaces. The complex phase of the eigenchannel, 
idicated and shown in color, is similar that of a Bloch 
wave at the Fermi level in the infinite atomic gold chain 
(half- filled 6s-band). 

Gold atomic wires has been realized and stud- 
ied experimentally by several different techniques. 
The low bias elastic and inelastic (phonon scatter- 
ing) transport is well characterized and understood. 



see Refi^ and references therein. The first eigen- 
channel (from left) at the Fermi-energy is shown in 
Fig. [3] for a 7-atom gold chain. Not surprisingly, the 
majority of the transmission is carried by the first 
eigenchannel (Ti = 0.994) with only a small trans- 
mission for the other eigenchannels (< 10~^). In ad- 
dition, it is clear that the current through the wire 
is carried by the 6-s electrons forming a half-filled 
one-dimensional band where the sign of the wave- 
function change by a factor of i (right moving) along 
the wire. The only difference with the corrsponding 
eigenchannel from the right (not shown) is the phase 
factor which correcsponds to a left-moving wave. 




FIG. 4: (Color online) Left eigenchannel for an OPV 
molecule bound by thiols to the hollow sites on Au (111) 
surfaces. The colors correspond to the two different signs 
of the almost real- valued wavefunction. 

For molecular wires, the experimental and the- 
oretical understanding of electron transport is less 
well understood. The calculated transmission 
through the OPV molecule shown in Fig [His 0.037 
and 100% of the transmission is carried through the 
first eigenchannel. Since the wave is almost totally 
reflected, the imaginary part of the wave-function is 
too small to be seen in the figure. In the calculation, 
the thiol bonds to the hollow site on the Au (111) 
surface and clearly shows that the conjugation of the 
molecule continues through the sulfur atom and that 
there is significant coupling to the gold leads. 

To investigate an asymmetric case, we carried out 
calculations on an OPE molecule bound by a thiol to 
the left lead, and with a tunneling barrier (hydrogen 
termination) to the right hand lead, see Fig. [5] The 
calculational details are the same as for the OPE 
molecule in Refi^. Because of the tunneling bar- 
rier, the transmission, Ttot = 0.0026, is lower than 
for the OPV molecule, and the left and right eigen- 
channels are considerably different. This can easily 
be understood by the large reflection at the right 
junction. 

In the three examples described here, we find that 
the symmetry of the eigenchannels can be intuitively 
understood from the band-structure of the corre- 
sponding infinite wires. For the gold-wire, the 5d- 
band is below the Fermi-energy which is situated ap- 
proximately at half-filling of the 6s-band. The cor- 
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FIG. 5: (Color online) Left(a) and right(b) eigenchannel 
for an OPE molecule strongly bound by a thiol to the 
left surface and weakly interacting with the right lead. 

responding infinite wires (polymers) for the molec- 
ular wires have energy gaps at the Fermi-energy. 
The eigenchannels therefore show the exponentially 
decaying solutions of the 7r-electron state in the 
complex-bandstructure at the Fermi-energy21. 

IV. SUMMARY 

We have in this paper developed a method to cal- 
culate the scattering states corresponding to elas- 
tic eigenchannels. The method is summarized in 
Eqs. (p8)) - ((32)l where the eigenchannels are found 
from quantites normally available in transport calcu- 
lations using the NEGF technique. In addition, we 
show three brief examples of elastic scattering states 
calculated for molecular and atomic-wires connected 
to three dimensional contacts. The eigenchannels 
for these systems can be understood from the band- 
structure of the infinite wires providing an intuitive 
understanding. 

The eigenchannels are useful to interpret elastic 
electron transport through junctions. We belive that 
they will be especially useful to investigate the effect 
of the contacts between device and leads, e.g., bind- 
ing site of the thiol-bond on Au-surfaces. In addi- 
tion, the method gives a useful basis to understand 
the effects of phonon scattering on the conductance 
and their propensity rule o^^i^^ . 
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APPENDIX A: APPENDIX: 
ORTHOGONALITY OF SCATTERING 
STATES 

Viewing the Bloch states in the infinite, periodic 
leads, as a starting point for perturbation the- 
ory, we can obtain the totally reflected solutions \ui) 
of a semi-infinite lead. In this case the perturba- 
tion is the removal of the coupling between the pe- 
riodic cells at the surface. Furthermore, the totally 
reflected states may again be used as the starting 
point in a perturbation calculation to obtain the full 
scattering states |^'/.r)- In this case the perturbation 
is the device region and its coupling to the leads. We 
will here show that the perturbation expansion gives 



solutions that are orthogonal and normalized. To do 
this we focus on the perturbation expansion of l^*;) 
from \ui) and note that the same derivation may be 
used to obtain the 1?!;) from \ui). 

Starting with a set of orthogonal and energy- 
normalized eigenfunctions |M„(i?)) of the isolated 
leads {n e /, r) we generate the full scattering states 

l^niE)) = G{E)V\un{E)) + \u„{E)) , (Al) 

where V — Vl + Vr. The response given by the 
retarded Green's function only contains waves trav- 
eling outwards from the device region. To show that 
the solutions generated in this way are normalized 
we use the Lippmann-Schwinger equation 

!*„(£;)> = Go{E)V\'fn{E)) + \un{E)) , (A2) 

where the unperturbed Green's function is Go (E) ~ 
(E-HL-HR-Ho + iS)-^. Together with Eq. (|JT|l . 
we obtain 



{^.4E)\^n'iE')) = {u,,{E)\^,AE')) + {uniE)\V^GHE)\^„,{E')) (A3) 
- {ur.{E)\u„,{E')) + {u„{E)\GoiE')V \^n'iE')) + {un{E)\V^ GHE)\^rAE')) (A4) 

= 6„^„,S{E E') + {ME)\v \^AE')) (^rrirT^ + ^^^r—s) (^5) 

= S„,n'SiE ^ E') , (A6) 

I 

which shows that the final scattering states are orthogonal and normalized. 
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The flux vi is the band- velocity divided by the L unit- 
cell length, while the left lead fc-vector is in units of 
the inverse unit-cell lengths. Likewise for R. Alterna- 
tively, k has a unit of inverse length while the v's are 
in units of velocity. In any case the unit of the flux 
normalized wavefunctions is square root of time. 
Since we are interested in scattering states around 
the Fermi-energy, we may choose an energy arbitrarily 
close to the Fermi-energy where there are no localized 
states. 



